Disorder-induced rounding of the phase transition in the large g-state Potts model 
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The phase transition in the g-state Potts model with homogeneous ferromagnetic couplings is 
strongly first order for large q, while is rounded in the presence of quenched disorder. Ifere we 
study this phenomenon on different two-dimensional lattices by using the fact that the partition 
function of the model is dominated by a single diagram of the high-temperature expansion, which 
is calculated by an efficient combinatorial optimization algorithm. For a given finite sample with 
discrete randomness the free energy is a piecewise linear function of the temperature, which is 
rounded after averaging, however the discontinuity of the internal energy at the transition point 
(i.e. the latent heat) stays finite even in the thermodynamic limit. For a continuous disorder, 
instead, the latent heat vanishes. At the phase transition point the dominant diagram percolates 
and the total magnetic moment is related to the size of the percolating cluster. Its fractal dimension 
is found df = (5-|--\/5)/4 and it is independent of the type of the lattice and the form of disorder. We 
argue that the critical behavior is exclusively determined by disorder and the corresponding fixed 
point is the isotropic version of the so called infinite randomness fixed point, which is realized in 
random quantum spin chains. From this mapping we conjecture the values of the critical exponents 
as /? = 2 - df, Ps = 1/2 and u = 1. 
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I. INTRODUCTION 

Quenched disorder has generally a softening effect on 
the singularities of phase transitions in homogeneous 
systems. As an example we mention that in the two- 
dimensional (2d) Ising model the ordered phase (and thus 
the phase transition) is washed out by random fields 
whereas due to layered bond randomness the Onsager 
logarithmic singularity in the specific heat of the ho- 
mogeneous system will turn into a very weak essential 
singularity (McCoy- Wu model) Singularities at first- 
order phase transitions are also soften^. In 3d the latent 
heat is either reduced by a finite extent (for weak dis- 
order) or the latent heat completely disappears and the 
transition becomes second order (for sufficiently strong 
disorder)^. In 2d, according to the rigorous result by 
Aizcnman and Wehii there is no phase coexistence in the 
presence of continuous disorder, thus the phase transition 
is always continuous. Thus there exists a set of random 
fixed points whose possible classification is a great chal- 
lenge of statistical physics. However, as in the theory of 
disordered systems, exact results in this field are scarce 
and therefore our present understanding of the subject is 
limited. Since a weak-disorder perturbation calculation 
starting from the pure system's fixed point is not feasible 
most of the existing results are numerical. In this respect 
an important role is played by the g-state ferromagnetic 
Potts model^, in which the homogeneous phase transi- 
tion is of first order for g > 4 in 2d'' and for (7 > 3 in 
3d. 

According to Monte Carlc^ and transfer matrix^ in- 
vestigations the critical behavior of the 2d random bond 
Potts model (RBPM) is governed by a line of g-dependent 



fixed points. For a given q the critical exponents are 
expected to be disorder independent, however gener- 
ally there is a strong cross-over regimeifi. The anoma- 
lous dimension of the magnetization, x, shows a smooth, 
monotonously increasing g-dependence, which saturates 
at a given finite value as g — *■ 00. On the other hand the 
critical exponent of the correlation length, v, which is 
related to energy-density correlations, is found to show 
only a very weak g-dependence. The measured values 
are all around the rigorous bound^^, v = 2/d = 1, which 
is required by general theorem for disordered systems. 
(Note however the confiicting result in Ref. 0, which is 
probably due to the use of the bimodal, i.e. discrete form 
of disorder.) 

Despite of the intensive numerical studies performed so 
far several aspects of the phase transition in the RBPM 
are not clarified yet. Here we mention that it is still little 
known about the mechanism which leads to the softening 
of the first order transition. The effect of different types 
of disorder, continuous or discrete, has not been investi- 
gated. Also universality of the transition with respect to 
the lattice structure is still an open question. Finally, and 
most importantly, there is no example for a specific sys- 
tem in which the critical properties are (conjecturedly) 
exactly known. 

In this paper we are going to investigate the aforemen- 
tioned questions using the example of the 2d RBPM in 
the large-(7 limit. The use of this model has several advan- 
tages, i) The homogeneous model, which has a strongly 
first-order phase transition, can be solved analytically-, 
ii) In the presence of randomness the model shows the 
generic features of disorder-rounded first-order phase 
transitions, iii) It allows for an important technical sim- 
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plification since the high-temperature expansion of the 
modeii^ (with arbitrary ferromagnetic couphngs) is dom- 
inated by a single diagram, whose calculation is reduced 
to an optimization problem^'', iv) Finally, this optimiza- 
tion problem is polynomial and there exists a powerful 
combinatorial optimization algorithm.^^ which works in 
strongly polynomial time, i.e. the necessary computa- 
tional power does not depend on the actual form of the 
disorder. 

This model, for which we shall refer simply as random 
bond Potts model in the following, has already been con- 
sidered in several previous work. Cardy and Jacobsen— 
has shown a mapping with the random-field Ising model 
in the level of interface Hamiltonians, which explains the 
vanishing of the latent heat at the transition point. Nu- 
merical studies for large (but finite) g-values^- predicted 
a smooth and non-singular behavior of the critical expo- 
nents as q —^ oo. Relation of the random bond Potts 
model with an optimization problem was introduced and 
studied by the method of simulated annealing in Ref. 0. 
The optimization problem was mathematically analyzed 
in Ref. Il4i Here its relation with sub-modular function 
optimization was shown and an efficient computational 
algorithm was developed. Some numerical results ob- 
tained by the optimization method has been briefly an- 
nounced in a LetteijM. 

The structure of the paper is the following. The model, 
its random cluster representation and the way how its 
thermodynamical and correlational properties are related 
to the dominant diagram are given in Sec.^ In Sec. IIIII 
the mechanism of breaking of phase coexistence is ana- 
lyzed and an estimate of the breaking-up length is given 
through extreme value statistics. Thermal quantities (in- 
ternal energy, specific heat) are calculated in Sec. IIVI In 
particular we study the effect of discrete and continuous 
disorder as well as the location of the transition point for 
different 2d lattices. The magnetization properties of the 
model, which are related to the percolative properties of 
the largest connected cluster of the dominant diagram is 
studied in Sec. El In Sec. IVII we show arguments and 
approximate mappings between our model and random 
quantum spin chains which leads us to conjecture the ex- 
act values of the critical exponents. Our conclusions are 
drawn in the final section while two simple examples are 
given in the Appendix. 

II. FORMALISM: THE LARGE Q-STATE 
POTTS MODEL 

The g-state Potts models is defined by the Hamilto- 
nian: 

n = -J2x,,6{a,,aj) (1) 

in terms of the Potts-spin variables, di = 0, 1, • • • , g — 1, 
at site i. The summation runs over all edges of a lattice 
S -E, and in our study the couplings, Jy > 0, are 



i.i.d. random variables. We consider the system in the 
large-g limit where, being the entropy per site s ~ Ing, 
it is convenient to introduce the reduced temperature, 
T' = T In g = 0(1) and its inverse as /?' = I/t' = /?/ In q. 
In the following we use these reduced variables but omit 
in the notation the prime. In terms of these variables the 
partition function is given by: 

Z = ^g-'3«(W) (2) 

which is convenient to express in the random cluster 
representation^^ as: 

GCE ijeG 

Here the sum runs over all subset of bonds, GCE and 
c(G) stands for the number of connected components of 
G. 

In the large-g limit, where ql^-^'^ ^ 1, the partition 
function can be written as 

Z = ^ g^(«) , <^(G) = c(G) + /3 ^ J,, (4) 

GCE ijeG 

which is dominated by the largest term, (f>* — maxg 4>{G), 
so that 

Z = no(7^*(l + ...). (5) 

Here the degeneracy of the optimal set G* is tiq = 0(1) 
and the omitted corrections in the r.h.s. go to zero for 
large q. The free-energy per site is proportional to 4>* 
and given by: 



where N stands for the number of sites of the lattice. In 
the Appendix we give an illustration of the convergence 
with q of the family of function fq (/?) . 

All information about the random bond Potts model is 
contained in the optimal set G* . The thermal properties 
are calculated from <j)* , whereas the magnetization and 
the correlation function are obtained from its geomet- 
rical structure. First we note that correlation between 
two individual sites is unity only if both are in the same 
connected cluster, otherwise the correlation is zero. In 
this way the average correlation function, C(r), is related 
to the distribution of clusters in the optimal set. Next 
we note that if, and only if, there is an infinite cluster 
the correlation function approaches a finite asymptotic 
value, limr_>ooG(r) = w?. The magnetization m, de- 
fined in this way, is the fraction of sites in the infinite 
cluster. In the ferromagnetic phase, T < Tc, the mag- 
netization is TO > and vanishes at the critical point, 
Tc, as TO ^ {—t)^, where the relative distance from the 
phase transition point is measured hy t = {T — T^jT^- 
As usual the divergence of the correlation length ^ at 
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the critical point is given in the form: ^ ~ \t\~'^, where 
^ is defined as the average size of the clusters. At the 
critical point the largest cluster is a fractal, its mass M 
scales with the linear size of the system L as M ^ U^^ , 
where d/ is the fractal dimension. Spin-spin correlations 
at the critical point decay algebraically as C(r) ^ r~^^, 
where the magnetization scaling dimension is given by: 
X = d— df = 13/v. 

In a system with a free surface one can define surface 
quantities, whose singularities are governed by surface 
exponents. In the ferromagnetic phase the fraction of 
surface sites belonging to the largest cluster defines the 
surface magnetization, which scales as: rus ~ (— 
At the critical point surface-surface correlations decay as 



by self- duality asi^: 



Cs{r) 



and the surface scaling dimension Xs is 



related to the surface fractal dimension dy: as d — 1 = 
dy- -\- Xs- We have also /3s = XgV. 

In the random bond Potts model in the large-g limit 
the structure and the fractal properties of the optimal 
set are dominated by disorder effects. Other examples of 
systems in which the disorder is dominant at the critical 
point are random quantum spin chaingli, for which the 
set of critical exponents is exactly known. This analogy 
and a possible isomorphism between the fixed points of 
the two problems, which will be discussed in detail in 
Sec. I VII bring us to conjecture the value of the critical 
exponents for the random bond Potts model, which we 
list below. 

The fractal dimension of the infinite cluster at the 
transition point and the scaling dimension of the (bulk) 
magnetization are related to the golden-mean ratio, (j) = 
(1 + v^)/2, as: 

df = l + 0/2, a; = 1 - 0/2 (7) 

The analogous surface quantities are conjectured to be: 

d} = 1/2, Xs = 1/2 (8) 

Finally, the correlation length exponent is given by: 



(9) 



thus saturating the rigorous bound for disordered 
systems >ii v >1. 

We close this section by a few technical remarks. The 
cost- function of the problem (j){G), given in Eq. (01, is a 
supermodular function, thus there exists a combinatorial 
optimization methodic to maximize it. For the particu- 
lar expression of 0(G) in the present model a specific ef- 
ficient algorithm, the Optimal Cooperation AlgorithmjAi 
has been formulated. The largest random system we 
treated by this method is a 512 x 512 square lattice. 

The investigations in this paper are restricted to two 
dimensions. In particular we work on the square lattice 
for which the location of the critical point is known for 
some distribution of the coupHngs. For < (3Jij < 1 
and for a symmetric distribution P{J + J) — P{J — J), 
where J is the mean coupling, the critical point is given 



T, = - = 27 



(10) 



In the numerical calculations we used either a discrete 
(bimodal) distribution: 

Pb{J) = p6(j - A - J) + (1 - p)S(j + A - J) (11) 

which is symmetric for p = I — p = 1/2, or a continuous 
(uniform) distribution: 



PuiJ) 



1/2A, forJ- 
0, otherwise 



A < J < J + A 



(12) 



In practical applications the latter distribution is ap- 
proximated by a set of discrete peaks (see for example 
Eq. (Iiyfl ). which is necessary since the algorithm works 
in terms of rational numbers. In the numerical calcu- 
lations we collected data about systems up to a linear 
size L = 256, The number of realizations were about 10^ 
for the smaller systems and around 10'^ for the largest 
one. To check universality, we considered non-self-dual 
randomness, analyzing triangular and hexagonal lattices 
with number of sites up to 16384. 



III. DESTRUCTION OF PHASE COEXISTENCE 
DUE TO DISORDER 

In this section we consider the phase transition in the 
non-random system and see how introduction of disor- 
der will destroy the phase coexistence, i.e. soften the 
first order phase transition into a second-order one. In 
particular in 2d we will estimate the disorder dependent 
length-scale 4, at which size the phase coexistence no 
longer exists. 

The solution of the optimization problem in Eq. 
depends on the temperature: the optimal set is gener- 
ally more interconnected in low temperature and con- 
tains more disconnected parts in higher temperature. In 
the homogeneous, non-random model with = J there 
are only two (trivial), homogeneous optimal sets: the 
fully connected diagram and the empty diagram. Con- 
sequently the free energy of the homogeneous system 
(/horn) is given by: 



iV/3/hom = 



H 



N(3Jz, 
T>T, 



T<T, 



(13) 



where z is the connectivity of the lattice. At the phase- 
transition point, f3c = l/Tc = (1 — l/N)/{Jz), there is 
phase coexistence, thus the phase transition is of first 
order and the latent heat, Ae, has its possible maximal 
value: f3cAe — 1. Note that this result holds for any 
regular (even finite) lattice. 

As quenched disorder is switched on, say in the para- 
magnetic phase, its fiuctuations will locally prefer the 
existence of the non-stable fully connected homogeneous 
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diagram, and this tendency is stronger around the phase- 
transition point, where the free-energy difference between 
the two homogeneous phases is small. This mechanism, 
as the temperature is lowered, will result in a sequence 
of non-homogeneous optimal sets with increasing fraction 
of connected component, whose average size is just the 
correlation length ^. If ^ stays finite at the phase tran- 
sition point then there is phase-coexistence and thus the 
first-order nature of the transition persists. This is the 
case in 3d for sufficiently weak disorder. 

In 2d, however, for any small amount of continuous 
disorder the correlation length is divergent at the tran- 
sition point thus the phase transition soften to second 
order. This is true in the thermodynamic limit, while 
in a finite system of linear size L weak disorder fluctua- 
tions could be not sufficient to break phase coexistence. 
The finite length-scale, Zf, = £, at which the breaking of 
phase-coexistence takes place can be estimated through 
extreme value statistics in the following wayS^. 

At the transition point of the homogeneous system, i.e. 
at /3 = l/(Jz), we compare the stability of the empty 
graph with a non-homogeneous diagram Gs, in which all 
sites are isolated except of a domain of S spins, which 
are fully connected. The relative cost-function is given 
by: 



ijeGs 

^ -^iC-Ja/A) + l (14) 

where we used that the distribution is symmetric, with 
J = J. Here in the second line the number of occu- 
pied edges of Gs is Sz minus the missing bonds at the 
perimeter of the fully connected domain, which is given 
by aS*^/^. In the last equation we used that we are at 
the critical temperature and introduced the representa- 
tion J2iji'^ij ^ J) = AS*^/^^, in terms of the Gaussian 
random variable C,, which has zero mean and variance 
one. 

Next we consider all independent positions of the fully 
connected domain in the finite lattice, whose number is 
denoted by Ns , and look for the largest value of the cost- 
function. Wo, in these diagrams. According to extreme 
value statistics^i in the large Ns limit its value is related 
to the cumulative probability distribution, P{W > Wq), 
and follows from the equation: 



NsPiW >Wo)^l 



(15) 



Making use of the relation for a Gaussian variable: 
P(C > Co) ^ exp(-Co/2)/Co, we obtain from Eq. ^ 
that Wo > 0, i.e. the new diagram has a larger cost- 
function, than the empty set, if Ns > exp(const(J/A)^). 
Here we are interested in the case when S is in the same 
order as L^, i.e. when there is a complete breaking of 



the homogeneous diagram, in which case Ns scales as 
Ns ^ i", with < a < 2. Thus the breaking-up length- 
scale is just lb ^ L and given by: 



lb « lo exp 



A 



(16) 



where Iq is a reference length. 




FIG. 1: Non-homogeneous optimal diagram for the bimodal 
distribution at a size L=256, corresponding to the breaking 
up disorder, A/J = 693/2693 = 0.2573. Tlie connected part 
of the diagram is compact and percolates in one direction. 

The relation between lb and A has been studied numer- 
ically, by calculating the smallest A at which the break- 
ing of phase coexistence in a system of size L takes place. 
(We used periodic boundary conditions.) An example of 
a non- homogeneous optimal diagram obtained in this way 
for the bimodal disorder is illustrated in Fig. 

For the bimodal disorder the distribution of the lim- 
iting values of A for a given size i = 64 is shown in 
Fig. 12 The distribution consists of several individual 
peaks, some of those are related to instabilities due to 
small finite plaquettes. The size of the smallest such pla- 
quette, I, is given by I ~ J/A, but the typical size of a 
connected domain varies between I and L, thus there are 
strong sample to sample fiuctuations. We have checked 
the validity of the relation in Eq. H16|) for the bimodal dis- 
tribution and the result is plotted in Fig.O In the region 
of size we could work L < 128 the asymptotic behavior 
in Eq. H16I) has not yet been reached, the data points 
deviate from the expected straight line. The effective be- 
havior of the breaking up length for these sizes is better 
described (the relative error is 0.003) by a dependence 
InL ~ J/A, as seen in the inset. In the numerical cal- 
culation J/A goes to infinity faster than the prediction 
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FIG. 2: Distribution of the breaking up strength of bimodal 
disorder for a size L — 64. The large individual peak at 
A/ J — 1/3 corresponds to an instability due to a 3 x 3 pla- 
quette. 



given above in Eq. because we implicitly consider 

in Eq. H14I) an euclidean domain S whereas the Optimal 
Cooperation algorithm finds more complicated domains. 
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FIG. 4: Distribution of the breaking up strength of quasi- 
continuous disorder composed of m = 127 individual peaks in 
Eq. 1191 for two different sizes. For larger sizes the separated 
peaks are due to the discreteness of the disorder. 
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FIG. 3; System size, L, as a function of the inverse average 
breaking-up disorder, (J/A)^, for bimodal distribution. (In 
the average of J/A we considereded the contribution of all 
peaks). For the sizes we considered a fit with (J/A) is better, 
as shown in the inset. 

We have repeated this calculation by using a (quasi) 
continuous disorder, which is approximated by 127 dis- 
crete peaks. As seen in Fig.0|now the distribution of the 
breaking strength A, for a relatively small size, L = 16, 
is Gaussian-like, however for a sizes L > 32, there is a 
secondary structure reflecting the effect of discreteness in 
the representation of disorder. As can be seen in Table I, 
the measured breaking strengths are consistent with the 
result in Ea. lfTC|l . unless for the two larger sizes, where 
the effect of discreteness of the distribution becomes sig- 
nificant. 

In the numerical calculations on the critical properties 
of the RBPM, whose results are reported in the following 



L 


<^> 






samples 


16 


0.34854 


0.00016 


0.3368 


31324 


24 


0.3275 


0.0002 


0.3408 


13927 


32 


0.3142 


0.0002 


0.3422 


14786 


48 


0.3133 


0.0002 


0.3801 


15561 


64 


0.3075 


0.0004 


0.3932 


2878 



TABLE L Breaking strengths for the continuous disorder (for 
the fourth column Lo — 1 is used). 



sections, we always choose the disorder strong enough, so 
that the relation L ^ lb is satisfied. 

At this point we comment on a relation between the 
RBPM in the large-;? limit and the random-field Ising 
model (RFIM) at T = 0, which has already been observed 
by Cardy and Jacobsen^. The form of the cost-function 
in Eq. ((njl is identical to that of the RFIM, where bonds 
in the RBPM are equivalent to sites in the RFIM having 
a random field of J^ — J and the perimeter contribu- 
tion in Ea. ll4|l is related to the domain- wall energy in 
the RFIM. Thus f3J corresponds to the ferromagnetic 
coupling in the RFIM. This mapping on the level of the 
interface Hamiltonian is however restricted to a smooth, 
non-fractal interface. Therefore it applies in the same 
way in the stability analysis of the pure phases, result- 
ing in a similar breaking-up length-scale as in Eq. H16|l . 
However the (fractal) structure of the clusters in the two 
problems are different. As a consequence the bulk and 
surface fractal dimensions (and the related critical expo- 
nents) are different in the two problems. For the RFIM 
the percolating clusters have the same fractal properties 
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as ordinary (site) percolation, as has been demonstrated 
by numerical calculations^?.. 



IV. SINGULARITIES OF THERMAL 
QUANTITIES FOR DISCRETE AND 
CONTINUOUS DISORDER 

In this section we analyze the singularities of thermal 
quantities, in particular we consider the internal energy, 
e, and the specific heat, cy- We start to rewrite Eq. © 
for the free energy of a given realization of disorder as: 

F = ~c{G*)T- (17) 

where the optimal set, G*, changes discontinuously as 
temperature is varied. For sufficiently high temperature, 
P' z max(J) < 1, G* is the empty set, whereas for low 
temperatures, min(J) > 1, it is fully connected. In 
between, around the transition point, however, there are 
more frequent changes. In a finite system of size L, the 
temperature interval of stability of the optimal set is 
found numerically to be around ^ 1/L, for the bimodal 
distribution. As a consequence in a finite system the 
free energy is a piecewise linear function of the tempera- 
ture, whose slope (the entropy, c{G*)), is monotonously 
increasing with the temperature. From Eq. 1)17(1 the in- 
ternal energy is obtained by derivation, which for a given 
sample is: 

E = - J2 J^o (18) 

This is a piecewise constant function of the temperature, 
thus in a given sample, even for a finite system the inter- 
nal energy shows discontinuities. We illustrate this point 
in a simple example in the Appendix, where the compar- 
ison between the bimodal and the uniform distribution 
for a general simple triangle is given. In the given (fi- 
nite) sample the phase transition point separates the two 
qualitatively different regimes in which the largest clus- 
ter is not percolating (paramagnetic phase) , from that in 
which it is percolating (ferromagnetic phase) . This phase 
transition in the given sample is sharp and of first order. 
Later we shall use this definition to locate the transition 
point of a series of samples and study their distribution. 

In the random system, as we mentioned before, the 
average quantities, such as the average free energy and 
its derivatives are of physical relevance. The averaging 
procedure generally acts to smear out the discontinuities 
in the internal energy. In this respect the behavior of 
the averaged quantities is different for the discrete and 
the continuous distributions. (See also the example in 
the Appendix.) As illustrated in Fig. |31for the bimodal 
distribution several discontinuities remain in the average 
internal energy, whereas for the uniform distribution the 
average internal energy is continuous, as illustrated in 
Fig. For the discrete distribution the origin of the 
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FIG. 5: Internal energy with bimodal disorder (/3Ji = 1/6 
and 13 J2 ~ 5/6) as a function of the inverse temperature, /3. 
Note on the discontinuities which are present after averaging 
in the thermodynamic limit, c.f. at /? = 0.75. In the inset 
the discontinuous behavior at the transition point, (3 — 5/6, 
is illustrated (the lines from bottom to top are for increasing 
size from L = 32 to L — 256). 

discontinuities is the degeneracies of (j) at some temper- 
ature. For example consider a square lattice where the 
bonds follow the bimodal distribution in Eq. 1111) . In this 
lattice let us focus on an elementary square, whose bonds 
are all strong. The sites of this plaquette are either non 
connected (giving a contribution to / of AT) or are fully 
connected (with a contribution to / of T + 4(J -1- A)). 
These two situations are degenerate at T = 4/3(J-|- A). 
The average jump in the internal energy, Ae, at this tem- 
perature can be estimated by the approximation with 
independent plaquettes, where a fraction of 1/8 squares 
have the above degeneracy, each of which giving a local 
jump per site of 3T/4. Thus Ae w 3T/32 which agrees 
reasonably well with the measured jump in Fig. |S| 

The average internal energy displays also a finite jump 
at the transition point, as illustrated in the inset to Fig.|31 
For the bimodal distribution this jump comes from those 
configurations in which a corner site has a strong and a 
weak bond, so that the isolated corner and the connected 
corner has the same contribution. The fraction of these 
corners is finite at the transition point, which leads to a 
finite average jump. Note, however, that the positions 
of these degenerate corners are distributed overall in the 
sample, thus their contribution is not exclusively related 
to the appearance of a percolating cluster. 

The jumps in the average internal energy will disap- 
pear, if a continuous distribution of disorder is used. To 
illustrate it we used a distribution, which consists of 2m 
discrete peaks of equal weight: 

m — 1 

/£ = 

-I- (5(7+A + ((m-l)/2-/c)e- J)] (19) 
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FIG. 6: Internal energy for (quasi-)continuous disorder as 
given in Eq. 1191 . In the inset the latent heat is plotted as 
a function of the inverse number of the discrete peaks in the 
distribution. 
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FIG. 7: Scaling plot of the singular part of the internal energy 
in Eq. 1201 for bimodal disorder. In the plot there is one 
free parameter, for which we used the conjectured value in 
Eq. ©. 



SO that in the largc-m limit we approach a continuous 
distribution. In Eq. H19|) e is a small number, and in 
the simulation it has been used e — 2J10^^. As shown 
in the inset to Fig. the jump in the average internal 
energy at the transition point goes to zero as ~ I/'ti, 
which can be understood as follows. For a given corner 
site out of different coupling sets there are only m 
for which the degeneracy, as described for the bimodal 
distribution holds. Therefore the average jump should 
scale with 1/m, as observed. 

Going back to the internal energy for the bimodal dis- 
tribution we argue that the true singularity can be ob- 
served at two sides of the transition point, i.e. as t ^ O"*" 
or i — > 0~. In a finite system of length L, we expect a 
singularity in the form: 

eit,L) ^ e{0+,L) +^-"€{11^/"), t>0 (20) 

and similarly for t < 0. Here the specific heat exponent, 
a, satisfies the hyperscaling relation: dv ^ 2 ^ a. The 
scaling function e(i/), is expected to be an odd function, 
thus, for small argument, to behave as e(y) ~ y. The 
scaling prediction in Eq. (|20|l has been checked and the 
results are shown in Fig. [7| The scaling collapse of the 
data is indeed very good with the exponents: i/ — 1 and 
a = 0, as conjectured in Eq. From the collapse of the 
points we estimate v = I within a numerical precision of 
1%. 

As we mentioned below Eq. (|I8|I one can define a tran- 
sition temperature, Tc(L), in a given finite sample, which 
separates the regime in which there is a percolating clus- 
ter from that in which the largest cluster is not perco- 
lating. (We note that the percolating cluster at strong 
enough disorder, i.e. when Ih <^ L, is fractal, in con- 
trary to the cluster at the breaking-up strength of dis- 
order in Fig. ^ which is compact.) The distribution of 
this temperature for bimodal and continuous disorder is 
shown in Fig. |H1 for periodic boundary conditions. For 



large systems one expects that the distributions approach 
a unique limiting curve in terms of the scaled variable, 
T — ic/cj where = T^i^L) — is the distance from 
the critical temperature in the infinite system and the 
variance of the distribution scales as <j{L) ~ L^'^ ^ with 
some scaling exponent, uj. As seen in Fig.|Hlthis limiting 
distribution, P{t), has its maximum in r > and the 
integrated distribution is : Pgp = Jq °° P{T)dT > 1/2. 
As a matter of fact Pgp is the spanning probability given 
by the fraction of samples in which the largest cluster 
is percolating. Indeed a large (but finite) sample perco- 
lating at T > 0, is also percolating at T = Tc- Accord- 
ing to our numerical results for the bimodal distribution 
Psp — 0.63(1). The finite size corrections to the spanning 
probability are given by: Pgp — Psp{L) ~ At{L), where 
At{L) ^ L~'^ is the shift of the critical temperature. Here 
we note that in a pure system the shift exponent is gener- 
ally given by: e = l/i^, if the transition is of second order. 
If, however the transition is of first-order the critical ex- 
ponent is given by its discontinuity fixed point valu o^^i^^ , 
which is Vd = 1/d. This is in accordance with the exact 
result in the pure model, in which the shift in the tran- 
sition temperature is given by ~ 1/-/V, as given below 
Eq. H13() . From the available data we could not make 
an independent estimate for e, however the data points 
seem to be consistent with e — 1/v — 1 for bimodal ran- 
domness. The exponent lu could be determined also with 
relatively large errors giving an estimate lu — 0.9(1). 

We have also studied the location of the critical point 
on the triangular lattice, for which the coordination num- 
ber is z = 3 and on the hexagonal lattice, which has 
z = 3/2. In Fig. O we present the distribution of the 
finite-size critical temperatures for different values of the 
number of sites N for the triangular and hexagonal lat- 
tice, with symmetric and continuous disorder and for pe- 
riodic boundary conditions. Here we take as reference 
temperature (3r = 1/(2 J), which is the transition point 
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FIG. 8; Distribution of the finite-size (percolation) transition 
temperature for the bimodai (a) and the continuous (b) disor- 
der for different sizes of the system. Note that the difference 
from the exact transition temperature is in the same order as 
the variance. 

for the non-random system and expected to be a good ap- 
proximation for symmetric disorder. The spanning prob- 
ability at the reference temperature is Psp = 0.65(5) for 
both kind of lattices, which is close to that in the square 
lattice. 



V. MAGNETIZATION AND THE FRACTAL 
PROPERTIES OF THE PERCOLATING 
CLUSTER 

Order and correlations in the RBPM are related to 
the structure of clusters in the optimal set G* . As we ex- 
plained in Sec. m at the critical point the largest cluster 
is a fractal, its mass M scales with the linear size of the 
system L as M ~ L'^^ . The structure of a typical opti- 
mal set is illustrated in Fig. 1 of Ref. ilgi for the bimodai 
distribution. Note that in G* each connected component 
contains all the possible edges, since the inclusion of any 
coupling with Jy > will increase the cost function. The 
topology of clusters in ordinary bond percolating is evi- 
dently different, since they contain all the strong bonds 
but none of the weak ones. 

We checked numerically that the largest cluster is in- 
deed a fractal by measuring the number of points in the 
cluster, /i(r, L), which are at most at distance r from a 
reference point. Similarly we measured the number of 
points in the cluster between a distance r -I- 1 and r: 
s{r,L) = fi{r + 1,L) — fi{r,L). Here we performed av- 
eraging i) over the position of the reference point and 
ii) over the disorder realisations. Under a scaling trans- 
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FIG. 9: Distribution of the finite size transition (percolation) 
temperature for triangular and hexagonal lattices with con- 
tinuous distribution of disorder for different sizes. 
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FIG. 10: Average mass of the percolating cluster within a 
region of size r at the critical point, for the square and the 
hexagonal lattices. The conjectured asymptotic behavior is 
indicated by a straigth line. Inset: average mass of surface 
points of the percolating cluster in the square lattice. 

formation, r r/b, fi(r, L) and s{r,L) are expected to 
behave as: 

s{r,L)^h'^f-^s{r/b,L/b) (21) 

Now in the first equation taking b — r we obtain, 
li{r, L) — r'^f jl{L/r), thus in a log-log plot one should 
obtain a straight line for L > r. As seen in Fig. IIUI this 
relation is indeed satisfied both for the square and for the 
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FIG. 11: Scaling plot of the cumulative distribution of the 
mass of the clusters at the critical point, using the conjectured 
value of the fractal dimension in Eq. l(7|l. In the inset the 
fractal dimension is presented as calculated from the ratio 
of the masses of the percolating clusters in different finite 
systems (see text). 

hexagonal lattices, in both cases - within the error of the 
calculation - having the same asymptotic slopes, which 
are compatible with the conjectured result in Eq. Q. 
We also performed a multifractal analysis showing only 
a single fractal dimension. 

To have a more precise estimate of the fractal di- 
mension we calculated s(r, L) for finite systems of size 
L = 16, 32, 64, 128 and 256. Comparing results of two 
consecutive sites we obtained df from the second equa- 
tion of Eq. H21|l by setting b — 2. The effective fractal 
dimensions as plotted in the inset to Fig. have only 
a weak r/L and L dependence, and from the asymptotic 
behavior one can estimate df = 1.81(1) in good agree- 
ment with the conjectured result in Eq. Q. 

Next we studied the cumulative distribution of the 
mass of the clusters, R{M,L), which measures the frac- 
tion of clusters having at least a mass of M. According 
to scaling theorjiS^ it should asymptotically behave as: 

R{M, L) = M-^k{M/L'^f ) (22) 

with T = {2 — df)/df. Note that the scaling relation 
in Eq. I|22|) contains only one free parameter, df, which 
can be estimated from an optimal scaling collapse of the 
distributions calculated for different sizes. As shown in 
Fig. the scaling collapse is excellent if the conjectured 
fractal dimension in Eq. Q is used. By analyzing the 
accuracy of the collapse we obtained an estimate: df ~ 
1.810(5) in very good agreement with the conjectured 
value. 

We close this section by studying the surface magneti- 
zation of the model. Using open boundary conditions we 
measure the average mass of the surface sites belonging 
to the largest cluster, Mg. According to considerations in 
Sec. ^ in a finite system at the critical point this should 



asymptotically behave as: Ms ~ L f. According to the 
data presented in the inset to Fig. ^1 the surface fractal 
dimension is given by = 0.495(10), which is in good 
agreement with the conjectured value in Eq. 

VI. ANISOTROPIC RANDOMNESS AND 
RELATION WITH THE STRONG DISORDER 
FIXED POINT OF RANDOM QUANTUM 
CHAINS 

In this section the random bond Potts model is re- 
lated to another problem of statistical physics, which has 
partially exact results about its critical properties. The 
analogies and relations between the models will be used 
to conjecture the values of the critical exponents in the 
random bond Potts model, which are already announced 
in Eqs.O,® and ©. 

The problem we consider first is the Potts model with 
correlated bond disorder, in which the vertical bonds are 
constant denoted by Jy , whereas the horizontal couplings 
have the same value in a given column, i = 1, 2, . . . , L, de- 
noted by Ji. Here the Ji are i.i.d. random variables and 
in the following we choose Jy = J. This model, which 
for Ising spins is the well known McCoy- Wu model^, is 
translationally symmetric in the vertical direction. Con- 
sequently the optimal set has a strip-like structure, and 
the 2d diagram G, is uniquely characterized by its Id cut, 
denoted by g. The cost-function of the problem can be 
written as 4>{G) = L^^'>p{g), where the size of the system 
in the vertical (horizontal) direction is iy ^ L (L) and 

i^ig) = ci(g) + -^C2(5)+/3 J2 -h+PJiL-c^ig)) , (23) 

is the cost-function of the Id problem. Here ci{g) is 
the number of isolated points and C2{g) is the number 
of connected diagrams which has at least two sites. In 
the thermodynamic limit the second term in the r.h.s is 
negligible and redefining the temperature we arrive to an 
equivalent cost-function: 

V;(5) = ci(g) + /3^ J,, , (24) 

with (3 = (3/ (1 — /3J) and at the critical point (3 J = l.The 
optimization problem in Eq. 1241) will be the subject of a 
separate publication. 

Next we consider the extreme anisotropic^ (time- 
continuum or Hamiltonian) limit of the strip-random 
Potts model, where we let J\\/J and the transfer 
matrix of the model in the vertical direction is written 
as: T — exp(— r?i). Here r is the (infinitesimal) lat- 
tice spacing and Ti. is the Hamiltonian of a random Potts 
chain in the presence of a transverse field?i: 

h 

7^ = -^J.<5(s„s,+l)-^-^M^ (25) 

i i ^ k=l 
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where Mi|si) = |si + l, mod g). If the couplings in vertical 
lines vary randomly from line to line then the transverse 
fields, hi, are random variables, too. The fixed point of 
the 2d strip-random Potts model and that of the random 
quantum Potts chain with Hamiltonian in Eq. (|25|l are 
isomorph, thus they have the same set of critical expo- 
nents. The latter system has been studied in detail by 
a strong disorder renormalization group method, which 
was originally introduced by Ma, Dasgupta and Hii2i and 
used later by Fishesii and others22i22i21i. This method is 
assumed to give asymptotically exact results, in partic- 
ular the critical exponents for q > 2 are identicalSS and 
are given hy^: 

i^±=2, j;_L = 1 - (l)/2, xl = 1/2 , (26) 

which also hold for the 2d strip-random model. 

Having the critical behavior in the presence of corre- 
lated, anisotropic randomness we try to relate it to the 
original problem, in which the disorder is isotropic. To do 
so we imagine that at the fixed point of the anisotropic 
problem we let the couplings to be random also in the 
vertical direction. In this way translational symmetry in 
the vertical direction is broken and in the originally ho- 
mogeneous strips connected and disconnected parts will 
appear. At the critical point these two (creation and de- 
struction) processes are symmetric, therefore it is plau- 
sible to assume that the mass of the largest cluster stays 
invariant, thus: AI ^ Ly x L^~^-^ ^ L'^f , from which the 
value of the fractal dimension, df — 2 — x±, as announced 
in Eq. lO follows. Repeating this argument for a surface 
cluster we obtain the result in Eq. Another, and re- 
lated assumption is that the correlation volume, i.e. the 
surface of the largest cluster grows in the same way in 
the two problems, as the transition point is approached. 
This means, that V{t) - (t)£,±{t) ~ L\\ \t\-''^ - ^(i)^ ~ 
jtl"^"^, from which ly = i^±/2 = 1 follows, in accordance 
with the announced result in Eq. 

VII. CONCLUSION 

One of the most spectacular effect of quenched dis- 
order is the rounding of first order transitions. In 2d 
this phenomenon can not be studied in a perturbative 
basis around the pure systems fixed point. In this pa- 
per we have demonstrated that an almost complete un- 
derstanding of this phenomenon can be achieved in the 
other limiting case, when the effect of disorder is strong 
and dominates the critical behavior. We argued that 
the appropriate system for this purpose is the random 
bond ferromagnetic Potts model in the large-g limit, in 
which thermal fluctuations are strongly suppressed. In 
the random cluster representation the properties of the 
system are related to the dominant diagram of the high- 
temperature expansion. At the critical point this dia- 
gram is a self-similar fractal, giving a natural explana- 
tion of the diverging length-scale, and its bulk and surface 
fractal dimensions are related to the critical exponents of 



the transition. We have introduced a plausible mapping 
between our isotropic system and that of strictly corre- 
lated disorder, which is isomorph with random quantum 
spin chains for which the critical propcries are known. 
The critical properies conjectured in this way are then 
checked by large scale numerical calculations based on a 
very efficient combinatorial optimization algorithm. We 
have studied in detail the possible differences in the re- 
sults obtained for discrete or continuous disorder. For 
discrete randomness, such as the bimodal one, which is 
frequently used in numerical calculations the internal en- 
ergy displays discontinuities, also the latent heat is finite. 
This observation should make people cautious while an- 
alyzing data obtained by the use of discrete disorder. 

The results of this paper can be extended in several 
directions. Including in the distribution of the couplings 
the values pjij < and f3Jij > 1 we obtain a system 
in which percolation and random bond effects compete. 
This problem can also be studied by our optimization 
algorithm and the properties of the stable fixed point 
can be calculated by analyzing the fractal properties of 
the relevant diagrams. Another extension is to include 
negative couplings, thus a Potts spin-glass problem is ob- 
tained. This still can be formulated as an optimization 
problem, which is most probably np-complete. Consid- 
ering the ferromagnetic case one can try to perform a 
large-g expansion around the conjecturedly exact results. 
For the critical exponents the corrections are expected to 
go as 1/ Ing^^:. Finally a further possible study is to clar- 
ify the dynamical properties of the system and to show 
if there are some relations with the ultraslow dynamics 
present in random quantum spin chainsi^^. 
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ucation under grant No FKFP 87/2001, by the EC Cen- 
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APPENDIX A 

Here we present two examples to illustrate the discon- 
tinuous nature of the internal energy in the g — > oo limit. 

We start with a 3 x 3 plaquette of the square lattice in 
which the couplings are 1 or 5 with the same probability. 
The realization we investigate is shown in the inset of 
Fig. To illustrate the approach to (7 ^ cx) we have 
calculated for different values of q the internal energy 
as a function of the reduced temperature (T — *■ Ting), 
as introduced below Eq.(^). As seen in Fig. E| for any 
finite q the internal energy is continuous and becomes 
discontinuous only in the g — s- 00 limit. 

Our second example is to illustrate the effect of round- 
ing due to averaging over disorder. We consider here an 



11 



40 
35 
30 



^25 

0) 



B 15 



A square 3x3 lattice with J=1 or J=5 with probability 1/2 



10 




P 



10 



q=2 

q=10„ 

q=10- 

q=10. 
q=10 

q infinite ■ 



12 



14 



16 



when c < a + b. We use the notation |a, b, c| = {a, b, c} 

and a < 5 < c (i.e. the variables with the tilde are 
ordered in ascending order). 

For the uniform distribution of disorder in Ea. (|12f) with 
J = A = 1/2 one finds for the average internal energy 
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FIG. 12: Average internal energy as a function of /3 for vari- 
ous values of q. In the inset the analyzed plaquette is shown, 
where the thick lines represent strong bonds (J = 5) and the 
thin lines weak ones (J = 1). 



elementary triangle in which the edges carry the (pos- 
itive) weights: a, b and c (a = pja, etc.). By simple 
inspection one finds that, for a given triplet {a, 6, c}, the 
free energy is: 



1 + f5{a + b + c} if a + 5 + c < I 
3 a ^ <a + b + c 



when i 



b < c and 



l + /3(a + 6 + c) if c< 
-PFa,b,cm^ { 2 + /3? if 2 , ^ 

3 if i < S + 6 



/3 



< C 



which is a continuous function of the temperature. 

By contrast for the discrete bimodal dsitribution in 
Eq.((Tll) with 7=1/2 and A = l/VS, which has the 
same mean and variance as the uniform distribution the 
internal energy is given by: 



E{T) = < 
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This presents five discontinuities of width 0.079, 0.296, 
0.158, 0.671, and 0.296 at the temperatures 0.317, 0.423, 
0.789, 0.894, and 1.183, respectively. 
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